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Q ■ Abstract 

' The QCD gluon equation of motion is solved approximately by means of the Green function, 

ly-^ , This solution is used to reformulate the Lagrangian of QCD such that the gluon propagator 

appears directly in the interaction terms of the Lagrangian. The nature of the interactions is 
discussed. Their coordinate-space form is presented and analyzed in the static, non-relativistic 
case. 

O^; 1 QCD Lagrangian and Equations of Motion 

The QCD Lagrangian density is pj = c = 1) 

Cqcd = -\F(-^^'^Fl{^J+J2i^U''^'W,,-m,6,,)i;i (1-1) 

in 

en ' where 

(N ■ i^^J = A^:) - A^^) - f,,,Af ^W, (1-2) 

and 

{D^)ij = 5ijd^+\gstf^ A^^) or = Isxs + i 5s 4^^) t('^). (1-3) 

The QCD coupling constant is Qs, fabc are the structure constants, and 

where the As are the Gell-Mann matrices. As usual, repeated indices are summed over, with the 
colour indices i,j = 1,2,3, the flavour indices q = d,u, s,c,b,t, and a,b,c = 1,...,8 for the gluon 
fields. Gluon indices will be usually indicated by round brackets to avoid confusion with vector 
indices; eg. a\^^ or A^"-^f^. 



It is convenient to write 



where 

is the free gluon field tensor and 
is the non-Abelian part. 



FirJ = ^tl + (1-5) 

Tf,=d,A^^ -d^Af (1-6) 
Gf, = -g,f^,,AfA^;}^ (1-7) 
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Upon substituting (jl-3p into (jl-ip . the Lagrangian density (|l-ip can be written as 

^QCD = + Cjr + Ci^ + Ci^ + Ci^ (1-8) 

wliere 

c^ = Y.^^'^^^^^-'^<i)'^i (1-9) 

/:^ = -i^('^)^- (1-10) 

-C/, = -j^;^'A^) where jj^'^ = <7. E^^^ff ^'^^^ (l-^l) 

£,3 = = i 5s /abc (S^^'^) - 5.4«)) A^,) = 9sfa,cAl^ Al^d,A^\ (1-12) 

= -7 G^"^'^'^ = -1^2 ^^^^ ^^^^ ^(c) (1.13) 



4 A**^ 4 

For the quark equations of motion, we need 

dC d ( dC 







and using p-9p and (|l-lip we acquire the colour Dirac equations: 

(i7'^9^-m,)V^-5stS?4^"V^ = 0. 
Equations p-lSp can be written as 

(I?^),^.V^,-m,V^ = 0, 

or, in matrix notation, 

where = {D^).^ is defined in equation (jl-3p and = [ipgjipgjtjjg] (T stands for transpose). 
The equation of motion for the gluon fields is 



(1-14) 



(1-15) 



(1-16) 



(1-17) 



(1-18) 



where t^") 

where 

and 



(cf. ref. [2j). In component form, eq. (jl-18p is (cf. ref. [3]) 

f ix) = ^'^^ (x) + ('^) (x) + f/^ (x) , 
q 

= gs fate (^^0 d.^h + 2 ^J) 9.^(0 + 4) ^^^4'^) ' 

Jg2 — Qs Jabc Jbde ^(rf) ^(e)- 



(1-19) 

(1-20) 

(1-21) 

(1-22) 
(1-23) 



2 Reformulation 



For the study of inter-quark interactions (and subsequently for the study of the properties of mesons 
and baryons) it is convenient to use a (formal) solution of the gluon equations of motion ()l-19p to 
reformulate the Lagrangian, and thus the action, of QCD, so that the gluon propagator appears 
directly in the interaction terms. Such reformulation has been shown to be useful for the study of 
inter-particle forces in scalar theory with a nonlinear mediating field [H |5] . 

The formal solution of (I1-19P involves the use of the symmetric Green function of that equation, 
and this requires a choice of gauge. We shall use the Lorentz gauge, d^Af^^"'\x) = 0, whereupon 
the "glue" equation (jl-19p can be rewritten as an integral equation, 

^/^ ^(x) = J dV D{x - x') f^^Xx'), (2-1) 

where D[x — x') \s a, Green function defined by 

d^d''D{x-x') = 5^{x-x'). (2-2) 
The Green function can be written as 

D{x - x') = 1 ^ D{k) e-^-(--'), (2.3) 

where D{k) =V/k'^ is the momentum-space representation of the Green function. (Note: We use 
the notation that x, k stand for 4-vectors, i.e. x = (t, x) and k = {k^, k), etc.) 

We have not included the free-gluon solution of the homogeneous equation (|l-19p in (|2-ip , since 
free gluons do not arise and so free-gluon solutions will play no role in the present considerations. 

The expression ()2-ip is only a formal solution of ()l-19p (in the Lorentz gauge), because the 

components jg ^""^ + of the current j'^ ^""^ depend on A^^ ^'^^ . Unfortunately, it is not possible 
to obtain an explicit, closed- form solution of the non- hnear equation ([2^ (i.e. of ([^49]) ) for yl^(") 
in terms of the quark fields V'g (at least we do not know how to do so). Thus, one must resort to 
approximation methods. 

Equation (12- ip can be solved as an iterative series for A^^"''^ (cf. ref. [6j). The first order term 
in this sequence is just the expression (12-ip . but with replaced by J^^"^ only, that is 

4(;)(x) = / dV D{x-x')j^^''\x') = I dV D{x-x') gsi>\{x')Yt^^^l,{{x'). (2-4) 

J q 

The index (1) in (l2-4p . as in all that follows, indicates that it is a first-order iterative expression. 
We note that the gluon fields are expressed here explicitly in terms of the quark fields V'g and the 
propagator D only. 

Correspondingly, to first order, this modifies the Lagrangian density (|l-8p . to 

LI = - + C,, (4f ) + (4f ) , (2-5) 

where, in (|2-5p . A^f^^^'^ is as given in (j2-4p . 

In obtaining LI, eq. (j2-5p . we have used the fact that, in the Lorentz gauge, 

Cr = -\d^A^^^ d^A^ ('^^ ~ \a^^^ d^d^A^ , (2-6) 

where ~ means equality modulo surface terms. Then, because dfj_d^ A'^ ^""^ = J^^"^ in first order, this 
means that Cjr = —A^'^j'^^^'^ (in first order). Thus, in light of ()2-4p . the reformulated Lagrangian 
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density LI (eq. (j2-5p '). and so the corresponding Hamiltonian and the action, is a functional of the 
gluon Green function D{x — x') and the quark fields ■0* only. We shall refer to the theory based 
the Lagrangian density LI as the "reduced model". 

At this point we might mention that in the QED (or SU{\)) case, the non-Abelian terms Cj^ 
and Ci^ do not arise, and LI corresponds to the reduced QED Lagrangian. This reduced QED 
Lagrangian was used previously (in the quantized Hamiltonian formalism) to derive relativistic few- 
fermion equations and, from them, the relativistic energy spectra for all bound states of positronium 
[7], muonium [8] and negative positronium and muonium ions [9j. The results are exact to O(q^), 
including, for positronium, the field-theoretic virtual annihilation correction. 

3 Analysis of the inter-quark interactions in the reduced model 
including the static limit 

We denote the action by 5 = j d^x C{x) and shall consider, in turn, each of the interaction terms 
of the action in the reduced model. Thus, the term corresponding to £/^, eq. (jl-lip . is 

^ {-\f/^i-)^W,^''^)=-\ j d^-ciV4(")(x)I)(x-x')4'^i(x), (3-1) 

where the quark-field current j^^'^\x) is given in eq. (also in (|l-2ip ) . Since J^^"^ ^ 9s, we 

see that Sj^^ corresponds to an energy contribution of 0{g1). We also see that this interaction term 
describes the quark currents interacting via the gluon Green function D, which, in the quantised 
theory, would correspond to a one- gluon exchange inter-quark interaction. 

To understand the physical content of the interactions corresponding to S^i^ {k = 1,2,3), it is 
useful to consider the case of static sources in the non-relativistic limit, as was done for the scalar 
model [HE]. For the static case, j'^^"'\x) = j'^^^\x), hence the static version of eq. p-ip . is 

sil^ = -\ I dtdt' d^xd^x'j';^''\x)Dix-x')jl;lix'). (3-2) 



5« = / 



Since 



D{x - x') = -6 ({x - x'f) = — ^— \6{t -t'-\x-x'\) + 6{t -t' + \x- x'\)] (3-3) 
vr Svr I a; — a; I 

it follows that 

J ^t'D{x-x') = lj^Jdt' [6it-t'-\x-x'\) + 6{t-t' + \x-x'\)]=l-j-^^. (3-4) 
Thus, in the static case, Sj^\ eq. (I3-1|) . can be written as 

5« = /dt L« = - |dt where = i | d^x d'x' j^^^\x) (^-S) 

It is clear from p-5p that H^^^ corresponds to the two-point potential energy function 

0^ 1 1 

V{xr,X2) = ^ r = — , (3-6) 

47r \xi — X2\ 47r xi2 

which is the non-relativistic limit of the one-gluon exchange interaction (in coordinate represen- 
tation). It reflects the 0{g1) "Coulombic" contribution to the interquark potential. Note that it 
depends only on the distance, xi2, between the points xi and 0^2, as expected. 
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Similarly, the component of the action corresponding to Cj^, eq. (jl-12p . is, in first order, 

5};^ = j d^x c'-^^ix) = -gsfabcj d^x ylg^(x)^j;5^(x)a^Ag^(x). (3-7) 

This term corresponds to an 0{g^) contribution to the energy. The Green function D appears 
in degree 3 in this expression, which reflects the 3-gluon vertex interaction corresponding to the 
"degree 3 in A^f''' form of d'j^ 



In the static case, = / d^x' D{x — x')j^^j^{x') becomes 



,) = I d'x'D{x-x')jl;l{x') = I dVig(:rO J dt'D{x-x') 
where we have used (|3-4p . and thus 



dv 

Att \x — X 



TV (3-8) 



47r •^^'^^^ ' dx^ \x - x'\ 



d?x' ■{a)i^l\ Xk X^ 



where A; = 1, 2, 3. Therefore, in the static case, 



j dtHf^ = -gsfabcj d^x A^^^{x)Af^^\x)dkA^^^''{x). 



(3-10) 



Thus, by §^ and (l3^ . 



H 



(1) 
h 



9s 



(4^)3 

gs 



fabc / <i^X 



\x — Xi\ 



d^xi 



,3 3%{X2) 
d X2 



kE - X2\ 



d'x3 4t(^3) 



\x - ajsl 



/abc y d?xi d?X2 d^X3 j^''^(a;i) j^''^''(a;2)4''^''(a;3) ^^fc(a;i, a;2, 2:3), 



where 



and 



(47r)3 

Uk{xi,X2,x^) = y d 

U^^\xi,X2,X^) = - j 



X — Xn 



_d_ 

\x — xi\\x — X2\\x — x^\'^ dxn 



d^a 



\X — Xi\\X — X2\\X — X^\ 



fnr 

J \v\\v 



U^^\xi,X2,X3), 



X2l\\v - X3i\ 



(3-11) 



(3-12) 



(3-13) 



here x^n = —Xnm = Xm — Xn- Note that the integral in ()3-13p is divergent and so must be 
regularised (this is discussed in ref. [3). However the derivatives (I3-I2p are well-behaved (finite). 
We see that the interaction term ()3-lip corresponds to a three-point potential, 



Vk{xi,X2,X3) 



{An 



-3 Uk{xi 



,X2,X3) = —^-^U^'^\xuX2,X3), k= 1,2,3, (3-14) 



which is an 0{gg) first-iterative -order "correction" to the two-point Coulombic interaction V{xi,X2) 
given in eq. (|3-6p . 

Lastly, the component of the action corresponding to Cj^, eq. ()l-13p . is, in first-iterative-order 

5};) = J d'x Cj,{x) = -\glfabcfadej d^X ^(x) .ig ^(x) ^jjj '^(x) '^(x), (3-15) 
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which is of 0{g^) and of degree 4 in D, corresponding to a four gluon interaction vertex. In the 
static case, this term becomes 

Sf^ = --/sfa.cfaaef d'x ^[J) ^(x) ^(x) ^(Jj '^(x) ^ig ^(x) = - j dt (3-16) 

where, by (|3-8p . 

- -/Ja.cfadej X A''!^ ^{x) A''^^ ^{x) A^^^ {x) A^^'' {x) (3-17) 



^2 



with 



C/(4)(a;i,^2,^3,^4) = r. r. r = rr^ n T' (^'l^) 

J \X — Xi\\X — X2\\X — X^WX — Xi\ J \V\\V — X2l\\V — X3i\\V — Xii\ 

If we include the couphng constants (and related factors) the corresponding potential-energy func- 
tion is 

1 

Vl^{xi,X2,X3,X4) = - U'^^\xi,X2,X3,Xi). (3-19) 

This is an 0{g^) first-order-iterative four-point potential correction to the Coulombic two-point 
potential (l3-6p . 

Unfortunately, the integrals ()3-12p . ()3-13p and ()3-18p that define the three and four point poten- 
tials cannot, in general, be evaluated explicitly, that is they cannot be expressed in terms of common 
analytic functions (at least we do not know how to do so). Nevertheless various general properties 
of these first-iterative-order non-Abelian corrections to the Coulombic inter-quark potential can be 
readily established, and analytical expressions can be obtained for particular situations. 

The general properties and representations of C/^ and C/*^^^ are presented and discussed in the 
Appendix, where it is shows, inter alia, that 

1) . The "potentials" U^'^\xi,X2,X'i) and [/(^) {xi,X2,X'i,X/i) are, in fact, functions of the 
distances Xmn = \xm — only (as might be expected of a closed system). 

2) . The evaluation of Uk{xi,X2.,x-i) can be reduced to the computation of two single quadra- 
tures and that of C/*^^-*(ki, X2, x^, X4) to a double quadrature, which must be done numerically. 

As mentioned, the three and four point functions f/^^^ and f/^^^ can be evaluated analytically 
for some particular cases. Thus, when all three distances are equal, i.e. X12 = X13 = X23 = r, 
the (regularised) three-point function of eq. (j3-14p becomes U^^\r) = 47rln(r/a), where a is an 
arbitrary distance scale [11 15]. Similarly, when points X2 and x^ are coincident, i.e. X23 = and 
2^12 = xi3, U^^\xi2) = 47r ln(xi2/a), and so 



Ukixi,X2,X3 = X2) = dU^^\xi,X2,X3)/dx^ = \dU^^\xi,X2,X2)/dx^2 = -2vr 



[Xi - X2)k 

333=332 kci - X2r 



This shows that the corresponding correction to the Coulombic one-gluon exchange potential (|3-6p 
due to the "cubic" interaction term (j3-12p . in the non-relativistic limit, is of the form (cf. (j3-14p ) 

Vi(x,x..X3 = x.) = ^(|^. (3-20) 

Similarly, for the particular case xi = X3 and X2 = X/j^, there is only one distance, X12, between 
the two pairs of coincident points, and the four-point potential function U^'^\ eq. (|3-18p . and so 
Vj^, eq. p-19p . can be evaluated explicitly: 

u^^^(^^^)= f r^2r-v — ^ = — ' Vi,{xu) = ^—. (3-21) 

J \V\^\V + Xi2r Xi2 A^-K X12 
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Although the expressions (j3-20p and (j3-2ip are only segments of the three and four point potentials 
Vk{xi, X2, X3), eq'n (j3-14p . and Vi^{xi,X2,X3,X4), eq'n (j3-19p . they suggest that their behaviour 
is Coulomb- like in general. Note that these corrections are 0{gg) and O^g^) respectively. 



4 Concluding remarks 

We have used an approximate, iterative solution of the non-linear classical equations of motion of 
QCD to derive expressions for the interaction terms corresponding to the non-Abelian terms ()l-12p 
and (jl-13p of the QCD action. In first iterative order, cf. equation (|2-4p . these turn out to be 
expressions involving products of three and four one-gluon exchange Green functions, corresponding 
to three- and four-gluon interaction vertices (cf. eq'ns (j3-7p and (j3-15p respectively). 

We have examined these non-Abelian terms in the static, non-relativistic limit and found them 
to be three- and four-point static potentials, (|3-14p and (|3-19p . that depend on the inter point 
coordinates only. Although we could not express these three- and four-point potentials in terms 
of common analytic functions in general, we could do so for some restricted sections and these 
indicate that the potentials Vk{xi,X2,X3) and V^^\xi,X2,X3,X4) are Coulomb-like in general. 

The derived three- and four-point cluster corrections together with the one-gluon exchange 
interaction could be used as a short-range contribution in potential models of baryons, tetra-quarks 
etc. 

The quantized theory, based on the reduced Lagrangian (j2-5p . can be used to derive relativistic 
few-quark equations, as was done for the scalar theory with non-linear mediating fields [HIS]. This 
shall be left for future work. 

Appendix. Properties and evaluation of the three- and four-point 
potentials 

Components of the three-point potential ()3-12p form the vector potential: 

d 

U{xi,X2,Xs) = —U^^\xi,X2,X3), (A-1) 

0x3 

where the function U^^\xi,X2,X3) (see eq. (|3-13p ) is studied in [4]. Particularly useful is the 
following representation this function: 

U^^\xi,X2,xs) = -^/d'^t Jd'x ^-kl{^-^i?-kl{^--2?-kl{^--z? 



'd?k_ 
W 



l"-! '^2-''12^'^2'^Z'''2'i^'^l'^3-''lZ/ 1 



00 



-Jd^k J —— e (^1^2^12+^2^3^23"'"^1^3^13)'^^ ^ (A-2) 




where Xmn = \xm — Xn\ {m,n = 1,2,3), k = k/k, k = \k\ = ^/kl + k2 + k^ and / d?k denotes 
integration over the unit sphere in 3D fc-space. The integral (|A-2p (as well as (|3-13p ) is divergent and 
needs to be regularized. One way is to split it into two terms: C/*-^-* = U + Uq, where C/(xi2, X23, X13) 
is a finite function of three scalar arguments and Uq = U (a, 6, c) is an "infinite constant" (a, 6, c are 
arbitrary constants) [4]. 
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Inserting ()A-2p into r.h.s. of (|A-ip discards the infinite constant Uq and yields the formula: 



U 



0[/(3) 



d 



dk 



^(1212^2 :i,2t.2^2 :i.2i.2„2 \u2 



2 [ d^kkl (kfx^i + kix32) I dk Aje-^'^i^i-L+' Ofc^ = x^j^ + 2:32/2, (A-3) 



where 



d k k^ k^ 



L2jL2^2 I 1,'z u'z 'z I KZ KZ 'z 
"'1"'2-^12 ' "'2'*'3'''23 "I" '''l'^3'''13 



21,2^2 



U2<,2^2 



n = 1,2. 



(A-4) 



Next, we introduce angular variables on the unit sphere in 3D A;-space, so that ki 

sin'!?cos(^, k2 = sin sin 93, k^ = cost?. Then 



2tv 



27r 



/ 



dip cos (f J, 



I' 



I2 = I ^(p sin^y? J, 




(A-5) 



where 



J 



TT 

/ 



di? sint?cos^ "!? 



(xi2 sin cos 99 sin ip)"^ + cos2t?(xf3 cos^y? + X23 sin ip) 



^12^^ 



sin 2lp R 
1 — arctan 



R 



sin 2p 



and i? = v^(cos 2p) + + r/2. 



"13 -^23 
2 ' 



T] 



2 _ [(a^l3 + a:23)^ - xl2][xl2 - {^13 " 2:23)' 



(A-6) 



(A-7) 



^2 -^12 

Inserting (|A-6p into (lA-Sh and using the integration variable s = cos 2(/7 yields the quadratures: 



h 
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'12 



'12 



d^ 

ds 



1 + s 1 + s i? 

arctan ■ 



1 - s 



1 - s 1 - s i? 

arctan ■ 



1 + s 



Note that for the special case rE23 = 0, the expressions (jA-SP to ()A-8P yield the result (|3-20p . In 
the general case, that is, for arbitrary values of x\2^ X13 and X23, the integrals ()A-8P need to be 
evaluated numerically. 

The behavior of the three-point vector potential as a function of 0:3 is illustrated in figure 1. 
We note the following symmetry properties of the vector potential: 

• translational invariance: U(x\ + A, + A, 013 + A) = \J(x\^X2^xj^^ where A G M^; 

• rotational covariance: L'"(Ra;i, Ra;2, Ra;3) = VJJ(x\^X2^x-^^ where R G S0(3); 

• partial permutational invariance: U{x2^X\^x-i) = XJ(x\^X2^x-^\ 

• scaling transformation: J7(Aa;i, Aa;2, Aa;3) = A^^[/(a;i, 0:2, 2:3), where A G M+. 
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Figure 1: Two-dimensional section of the vector potential U{xi,X2,X3) as a function of 2:3 = r 
for fixed Xi and X2- The arrows indicate the direction of the vector field U{xi,X2;r). Note that 
it is invariant with respect rotations about the axis xi-X2- The dashed line between xi and X2 
corresponds to J7 = 0. 



These properties follow from the properties of the three-point scalar potential (I3-13P stated in fUlS]. 
We also note that the scaling transformation of the potential U reflects its Coulomb-like behaviour. 



The four-point scalar potential (I3-18P can be treated in similar manner: 

J \X-Xi\---\X- Xi\ TT^J J 

00 



^jd^k Jdke--^ =J (A-9) 





where X = kik2Xi2 ■^13 ~l~ kik^Xi^ -\- ^2^3X23 -|- k2k^X24 -\- /C3/C4X34, 

kn = kn/k (n = 1, ...,4), k = ^yk1 + ■ ■ ■ + k^, and / d^^ denotes an integration over a unit hyper- 
sphere in 4D /c-space. 

Next, we introduce angular variables {x,'d,^} on the unit hyper-sphere in 4D A;-space, so 
that ki = sin X sin 1? cos 99, k2 = sin x sin sine/?, k^ = sin x cost?, k^ = cosX; and / d^/c = 
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J f siiiT^di? J sin^ xdx ■ Then using the integration variables u = cosx, v = cos??, w = cos 2^3 




reduces the integral ()A-9P to the form: 



1 1 




-1 



1 

where I = J 






= {ix?2(l - w^){l - v^) + i[x?3 + xl^ + (x^g - xl^)w]v^] (1 - v^), 
= i[x?4 + xl^ + {xj^ - xl^)w]{l - v^) + xl^v^ - 



This double integral can be evaluated numerically. Evaluation of ()A-10p for the case xi = and 
= gives the same result as (j3-2ip . 

Finally we note, that the four-point potential ()A-10p possesses the same translational invariance 
and scaling transformation properties as the three-point one does. Besides, it is invariant under 
arbitrary rotation and permutation of its arguments: 

• rotational invariance: U^^\Rxi,Rx2,Rx^,Rxi) = U^^\xi , a;2, a^3) 3^4), where R G S0(3); 

• complete permutational invariance: 

U^'^\x2,Xi,X2„Xi) = U^^\xi,X^,X2,Xi) = ■■■ = U''^\xi, X2, X^.X^) . 
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